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ABSTRACT 

Modern  Naval  weapon  and  sensor  systems  are  strongly 
influenced  by  the  marine  environment.  Foremen-  among  the 
atmospheric  effects  is  ducting  of  slectromagnatic  energy  by 
refractive   layers    in    the    atmosphere.  To    assess   the   effect 

of  ducting  on  electro iragnetic  amissions,  the  Havy  developed 
the  Integrated  Refractive  Effects  Prediction  System  (USPS). 
Research  at  Naval  Postgraduate  School  (NPS)  hss  led  to 
development  of  a  state-of-the-art  model  which  car.  he  used  to 
predict  changes  to  the  refractive  profile  of  'he  lower 
atmosphere.  The   model      uses   radiosonde      data    and      surface 

meteorological  observations  to  predict  changes  ii  refractive 
conditions  and  lew  level  cloud/fog  formation  over  18  to  30 
hour   periods.  The    model   shows      some    skill      in    forecasting 

duct  regions  when  subsidence  rates  can  be  specified  '•o 
within  +/-.0015  m/s.  This  thesis  shows  the  applicability  of 
the  NPS  marine  atmospheric  mixed  layer  model  to  flee-:  tac- 
tics. Atmospheric  refractive  effects  on  specific  emitters 
can  be  predicted  when  model  predictions  are  used  in  conjunc- 
tion   with    IREPS. 
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I.        INTRODUCTION 

Over  the  centuries  military  leaders  have  learned  to 
appreciate  and  take  advantage  of  the  effects  of  the  environ- 
ment, and  when  they  failed  to  do  so,  have  met  with  disaster. 
A=  it  approached  the  islands  of  Japan  in  the  thirteenth  cen- 
tury, Ghengis  Khan's  Mongol  invasion  fleet  was  destroyed  by 
a  typhoon.  From  August  to  October  15B8  the  Spanish  Armada 
was  beset  by  storm  after  storm  which  resulted  in  the  sinking 
of  many  ships  and  failure  of  Spain's  attempted  invasion  of 
Britain.  On  D-Day,  June  6,  1944,  the  allied  invastion  force 
took  advantage  of  predicted  good  weather  between  storms  and 
nade  •'-he  amphibious  assault  across  the  beaches  of  Normandy. 
Tn  December  of  the  same  year  Admiral  Haisey's  Third  Flee- 
was  caught  preparing  to  refuel  by  an  undetected  typhoon, 
causing  28  ships  to  be  crippled,  156  airplanes  to  b<=  lest, 
and  the  destroyers  S FENCE,  MONAGHAN,  and  HULL  to  be  sunk 
(Nash    1976). 

Modern  naval  warfare  technology  is  providing  increas- 
inaiy  capable  and  complex,  but  environmentally  dependent 
weapons/sensor  systems.  Net  only  the  severe  sea  and  weather 
conditions  must  be  predicted  t3  retain  advantage  over  adver- 
saries,        but    also      environmental    factors      which    enhance      or 


degrade  a  wide  range  cf  weapon  and  sensor  systems  which  uti- 
lize electromagnetic  (EM)  propagation.  Electromagnetic  fre- 
quencies abovt-  the  H?  (3-30Mhz)  band  can  be  greatly  affected 
by    atmospheric   refraction. 

K  condition  known  as  ducting  occurs  when  a  refractive 
layers  oause  EM  energy  tc  bend  toward  the  earth  at  a  rate 
greater  than  or  equal  to  the  earth's  curvature.  Ducting 
occurs  with  certain  critical  vertical  gradients  of  tempera- 
ture and  humidity,  and  can  cause  both  increased  radar  and 
radio   ranges,    and   holes    (gaps)    in   normal    coverage. 

Tactical  advantages  exist  by  knowing  duct  locations. 
These  e.5  ve  nta:jes  include  being  able  to  make  realistic  esti- 
mates of  "SM  detection  and  counter  detection  ranges.  This 
would  enable  commanders  to  mafca  decisions  concerning  emis- 
sion control  (EMCON)  ,  the  positioning  of  both  air  and  sur- 
face surveillance  assets,  the  altitude  and  flight  profiles 
for  strike  aircraft  to  minimize  detection,  the  placement  of 
electronic  jammers  for  maximum  effect,  and  numerous  other 
tactical    considerations. 

Ducting  commonly  occurs  with  inversions  which  act  as 
trapping  layers,  refracting  or  bending  EM  energy  toward  the 
earth.         Inversions   are   stable   layers   between    warm,      dry    air 
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aloft  and  cooler,  more  moist  air  below  and  typically  exist 
in  marine  surface  high  pressure  regions.  Extensive  areas  of 
low  level  stratus  clouds  often  delineate  areas  of  duct 
occurren  ce . 

Ducting  is  expected  to  be  minimal  near  fronts  and  areas 
of  convective  cloudiness.  Fronts,  with  their  associated 
upward  motion,  often  dissipate  the  inversions  as  the  whole 
air    column      becomes    mixed.  Areas   of      convective    activity, 

detectible  by  the  presence  of  oumuius  clouds,  are  also  nor- 
mally inversion  free.  Thus  satellite  infrared  (IP)  and  vis- 
ible photography  should  provide  means  to  estimate  ducting 
regimes. 

The  Navy  employs  a  microcomputer  based  system,  IREPS 
(Integrated  Refractive  Effects  Prediction  System)  (Hitney, 
1979),  to  identify  ducting  conditions  and  to  assess  the 
effects  or.  various  fleet  EM  emitters.  IREPS  reguires  radio- 
sonde  data    as   input. 

The  disadvantage  of  IREPS  is  -hat  it  predicts  the  duct- 
ing conditions  only  at  the  radiosonde  launch  site  and  only 
at  the  sounding  time.  As  such,  it  is  not  an  IREPS  weakness 
but  more  an  inherent  weakness  in  single  station  assessments. 
The    variability   of   the  atmosphere   makes    projections    in    space 
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and  predictions  in  tins  of  the  ducting  environment  important 
for  the  tactician,  while  satellite  data  can  provide  some  of 
the  data  necessary  to  make  projections  in  spacp,  something 
else  is  needed  to  predict  changes  with  time.  Considering 
EMCON,  the  ability  tc  forecast  ducting  for  a  24  hour  period 
may  be  critical.  Onder  strict  EMCON  conditions  radiosondes 
cannot  be  launched  because  the  signal  they  emir  while  trans- 
mitting  data    could   act   like   a   beacon    to   hostile    fcrcss . 

A  predictive  model  exists  for  changes  in  the  marine 
atmospheric  mixed-layer.  The  nicro  computer  based  model  was 
developed  by  the  Environmental  Physics  Sroup  at  the  Naval 
Postgraduate  School  (Davidson  st  al,  1933).  The  aodel  is 
initialized  with  IREPS  sounding  data,  surface  observations 
of  sea  surface  temperature,  wind  speed,  and  subsidence  at 
the    inversion. 

Utilizing  the  HP-9845  (which  ail  carriers  currently  have 
on  board  for  use  as  tar-  of  the  IREPS  system)  and  data  in 
IREPS  format,  Brower  {198  7}  developed  a  program  to  incorpo- 
rate the  NPS  mixed  layer  model  to  predict  changes  in  the 
refractivity      profile    for      a      33       hour    period.  The      model 

reguires  as  inputs:  surface  wind,  sea  surface  temperature, 
and    an   estimated   subsidence   value      as    well    as      current    IREPS 
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sounding  data  to  predict  ths  mixed  layer  evolution  and 
resulting  changes  in  both  surface  based  and  elevated  duct 
haighxs    up    tc    1200    m    altitude. 

The  purpose  of  this  thesis  is  tc  svaluate  this  state  of 
the  art  atmospheric  mixed  layer  model  and  its  application  to 
fleet  tactics.  The  data  utilized  in  this  study  were  from 
radiosondes  taken  by  USS  NIMIIZ  and  OSS  COEAL  SEA  in  the 
Arabian  Sea  in  xhe  vicinity  of  2  3N  361E  ,  a  data-poor  area, 
during  February    1980. 
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Figure   1.      Regional  Chart    Showing   Sounding    Data   Locations. 


13 


II.        THE    MODEL     DESCRIPTION 

A.        MODEL    OVERVIEW 

The  NPS  mixed  layer  model  is  an  integrated,  two  layer, 
zero-order  model  which  predicts  changes  i.i  the  marine  atmos- 
pheric boundary  layer  (MABL)  bslow  the  inversion.  Ihe  two 
layers  consist  of  the  lower,  wall-mixed,  turbulent  boundary 
layer  topped  by  the  relatively  non-turbulent  free  atmos- 
phere. The  two  layers  are  5= para ted  by  an  inversion  or 
transition  zone  which  is  assumed  to  hive  zero  thickness* 
The    predicted   properties    are: 

(1)  height    of  inversion 

(2)  values   of  weil-mixad    properties 

(3)  values   of   jumps    at    the    inversion 

(4)  formation  of    clouds/fog    within    the    mixed    layer 
Procedural    uses      are    shown   in      Figure    2,       and      the    equations 
will   be    presented    in    Section    3.       Inputs    to    the    model    srs; 

(1)  Radiosonde    data 

(a)  Vertical   distribution    of    temperature 

(b)  Vertical   distribution    of    moisture 

(2)  Sea   surface    tempera tura 

(3)  Wind    speed 
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(4)  Subsidence 

(5)  Latitude 

(6)  Julian  date 

(7)  Local  time  of  sounding 

From  the  predicted  properties,  a  ref racti vit y  profile  is 
produced  using  *he  refractivity  equation  (Section  B) - 
Observed  and  predicted  sea  surface  temperature  and  wind 
speed  are  input  parameters.  Subsidence  rites  are  estimated 
from  previous  observations.  Latitude,  Julian  da^-e,  and 
local    times   are    used    to    estimate    incident    solar   radiation. 
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Figure  2.   NPS  Mixed  Layer  Medal  Functional  31ock  Diagram 
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B.        EASIC    MODEL    EQUATIONS 

Once  initial  inputs  have  been  made  and  the  vertical 
structure    digitized   by  fitting   straight   lines   to   approximate 

the  vertical  profiles  of  temperature  and  humidity,  Equations 
1A,  1B,  and  1C  are  used  to  calculate  tha  flux  scaling  param- 
eters. The  calculations  are  based  on  the  air-sea  tempera- 
ture and  specific  humidity  differences  and  the  wind  speed. 
The  scaling  parameters  are  then  used  to  estimate  surface 
fluxes  of  momentum,  heat,  and  moisture  in  Equations  2k,  2B, 
and    2C. 

Equation  3A  is  used  to  predict  the  change  in  mixed-layer 
properties  by  the  difference  between  surface  and  inversion 
fluxes   scaled    by   inversion   height.  Equation    3B    is    used    to 

predion  changes  in  the  inversion  height  by  adding  the  flux 
at  the  inversion  scaled  by  the  inversion  jump,  a  measure  of 
the  strength  of  the  inversion,  to  the  subsidence.  Equation 
3C  models  the  change  in  the  inversion  ■jump  by  combining  the 
effects  of  surface  and  inversion  fluxes  with  subsidence  and* 
the  change  cf  the  inversion  height  determined  in  equation 
3B.  This  requires  use  of  the  gradient  of  the  property  above 
the    inversion. 
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Flux    Scaling    Parameters: 

u*  ■  ca^u 


T*    =    (9   -    9o)     Ce 


q*    =    Ce'*  (q    -    qc) 


Surface    Fluxes: 


-  (0  »W')    =    U* 


-  ( H  •  T '  }    =    Q  *T  * 


(1A) 
(1B) 
(1C) 

(2A) 

(23) 
(2C) 


Where: 


-  ( W  »  q  ; )     =    a  *q  * 

Cd,Ce   =   stability    dependant   draq   coefficients 
0*    =    wir.d   speed 

U*,T*,q*    =   scalinq    parameters    of    vertical    turbulen- 
momentum,    temp,    and   moisture  transfer 


-  (U '  W  •  )    =    downward    turbuien-    transfer    of    momentum 

- ( W ' T ' )    =   upward   turbulent   transfer    of   heat 

- ( W ' q • )    =    upward    turbulent    transfer    of    moisture 


Predictive    Equations     (general    form): 


dXm/dt    =     (W'X'c  -    W'X'i    *    f  (R)  )     /    h 


<3&) 
(3B) 


-W'X'i    =    DX    (dh/dt    -   Ws)     +    f(R) 

d(DX)/dt    =    Yx(dh/dt-Ws)    -     (W^X1  o-WMC1  i  +  f  (R)  ) /h    (3C) 


Where: 

Xm    =   any    well-ndxed   property:     0,    9v,    or   q 

9    =    potential    temperature 

q    =    specific   humidity 

H'X'    =    turbulent    vertical    transfer     (flux)    of   X 

o    subscript   =    "at    the    surface" 

i    subscript   =     "at    the    inversion" 

h    =    mixed-layer  depth    (or    inversion    height) 

DX    =   change    in    X    across   the   inversion:     (Xabove    -    Xm) 

Hs    =   subsidence  at    top   of    inversion 

Yx    =   vertical    gradient    of    X   abovs    inversion 

f  (R)    =    radiation    factor 

Pefractivity   Equation: 

M  =    77.6    P/T    ♦    6    X     105      (q?/T*)     +    .157    %  (4) 

Where: 

M  =    modified   r efr activity 

P  =    atmospheric   pressure    in    millibars 

T  =   temperature  in    Kelvins 

Z  =    altitude    in  meters 


19 


1000 


Q/kg 


25 


30    350  400 
M  UNITS 


Figure   3 


r  :.q. 


Vertical   Structure   of    9,    g   and    Corresponding   M 
Profile   for  a    Typical    Inversion     (tfGdel 
Simplification    is   Dashed)  . 

I      illustrates    the     relationship    between      vertical 


distribution  of  temperature  and  humidify  and  the  modified 
refractivity,  ?..  M  profiles  are  used  as  model  output 
because  of  the  ease  with  which  ducting  information  can  be 
extracted  from  them.  The  top  of  a  duct  corresponds  to  the 
height  above  the  surface  where  the  M  value  is  a  minimum. 
The  duct  base  corresponds  to  the  height  at  which  a  vertical 
line  drawn  dcwnward  from  the  point  of  minimum  M  value  first 
intersects   a    point   of   equal   M   units   or    tha    surface. 


20 


C.        SATELLITE    DATA    INFUT    TO    MODEL 

Provision  is  made  in  the  mDd<=l  for  incorporating  fore- 
cast wind  and  sec.  surface  temperature  changes  during  the 
prediction    period.       Although   ths    mixed   layer    is    sensitive    to 


SST, 


believed   that   accuracies   of    «•/-    1°C     in    5ST   are 


sufficient  foi  reasonable  model  accuracy.  Stewart  (198  1) 
reports  the  following  satellite  measurement  capability  has 
been  achieved:  NOAA-6  (VHRR)  infrared  radiometer:  SST  to 
.6°C  (with  no  clouds)  ;  SeaSat  and  Nimbus-7  microwave  radiom- 
eter (SHME)  :  SST  to  1.0°C  (with  no  rain,  no  RFI  (radio  fre- 
quency  interference),    and    >600    km   from  land). 

Satellite  SST  could  be  usei  in  the  flax  scaling  equation 
(13)  to  determine  9c.  Three  methods  are  currently  used 
aboard  ship  for  SST  determination:  Bathythermographs,  sea 
water  injection  temperature,  and  the  ' bucket /thermometer ' 
method.  None  of  these  yield  the  actual  sea  surface  "skin" 
temperature  which  is  the  relevant  quantity  in  the  determina- 
tion of  stability  and  turbulent  heat  and  moisture  fluxes. 
Therefore,  satellite  measurements  would  not  only  provide 
cloud  descriptions  over  broad  areas,  but  would  provide  a 
means   to    obtain   the   "skin"   temperature. 
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Ill,        ARABIAN   SEA    CLIMATOLOGY    AND    SYNOPTIC 
ITESCRTPTTUfl 


A.       CLIMATOLOGY 

The  climatology  of  the  northern  Arabian  Sea  is  dominated 
by  the  summfer  tsou-:hwest)  and  winter  (northeast)  monsoons. 
The  data  used  in  this  study  were  from  February  during  the 
latter  parr  of  the  northeast  monsoon  whi-:h  is  characterized 
by  light  offshore  winds  averaging  5-10  fe  riots,  and  normal 
Hadiey  circulation  with  the  subtropical  jet  located  near  30 
N  at  the  2  00  rab  level.  Fig.  4  shows  the  average  per  cent  of 
time  that  elevated  ducts  occur  during  February  to  March  in 
this  region. 
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Figure  4 


70"  E 


Freguency    cf   Elevated    Duct   Occurrence    Durinq 
Spring,    after    Lammers    et   al    (1930). 


FNOC       (Fleet    Numerical     Ocea nographic    Center)         products 


ware  use 3  to  c 


escribe  the  synoptic  situation  during  the  da 


period,  They  were:  Surface  Analyses,  530  mb  analyses  with 
Temperature  fields,  300  mb  analyses,  250  mb  wind  fields. 
Also*  NOAA-6  satellite  photos  were  used  to  determine  cloud 
coverage  cf  the  area  as  well  as  to  oonfirn  jet  positions. 

The  synoptic  situation  in  the  region  from  6-23  February 
1980  was  characterized  by  a  stable  high  pressure  surface 
ridge  ever  the  Arabian  peninsula  and  Arabian  Sea,  with  low 
pressure  troughs   over  the   Red  Sea  to   the  west   and  Indian 
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subcontinent   to      the    east.  The    most      significant   synoptic 

change  occurred  from  11  to  14  February  as  an  853  mb  low 
moved  across  the  northern  part  of  the  region.  This  low  was 
northwest  of  the  Persian  Gulf  a t  39  degrees  North  latitude 
on  11  February  and  tracked  eastward,  arriving  over  northern 
India  on  14  February.  It  was  during  this  brief  period  of 
11-14  February  that  the  minimum  shipboard  surface  pressures 
for  the  entire  period  were  recorded.  The  minimum  recorded 
surface    pressure    was    1012.0   mb   on    1'4    February. 

Table  I  lists  cloud  coverage  over  the  immediate  area  of 
the  data.  Low  level  stratus  clouds  were  confirmed  only  on 
13    February  and    from    19    to   21    February. 


TABLE    I 
Data    Area   Cloud    Coverage    9-21    February    1980. 


DATE  CLCUD    TYPE    PRESENT 

8-11  clear    with   scattered    cirrus 

11-12  no   satellite    coverage 

13  low   level    s-ratus 

14-18  clear    with    scattered   cirrus 

19-20  cirrus   and    stratus 

20-21  cirrus,    stratus   with   multilevel   to    NW 


^determined    from    NOAA   5    IR   imagery 


?ti 


Jin  a  tt  «3mp*  wis  made  to  derive  subsidence  values  from 
synoptic  scale  wind  fields  by  calculating  divergence  at  var- 
ious levels.  It  rfas  noted  that  at  upper  levels,  wind  maxima 
or  jet  streams  had  a  dramatic  influence  or  divergence  pro- 
files. Although  the  magnitudes  of  vertical  motion  obtained 
in  this  Banner  w  =  re  inconsistent  and  highly  variable,  an 
interesting  and    possibly    significant    observation    was    made. 

The  subtropical  jet  associated  with  the  wind  maxima  dis- 
cussed above  is  usually  found  at  the  northern  limb  of  the 
Hadley  Cell  (Fig.  5).  It  is  at  this  northern  limb  where  the 
greatest  downward  motion  or  subsidence  is  found  (Palmen  and 
Newton,      1969).  Although  this    is      a    feature   of      the    large 

scale  circulation,  it  seems  evident  ^hat  the  jet  does  have 
an  effect  wvich  should  be  considered  in  single  station 
assessments   ^hen^ver    possible. 

Changes  in  the  subsidence  rates  caicilated  by  hindcast- 
ing  and  shown  in  Table  II  correlated  closely  with  changes  in 
the  position  and  strength  of  the  upper  i  =  vel  jet  stream. 
Daring  the  period  6-9  February  there  was  a  steady  decrease 
in  subsidence  rates  from  -.0073  to  -.0355  m/s.  During  the 
same  period  the  tail  of  the  250  mb  jet  max  moved  steadily 
eastward    from    a    point   due    north    of   the    data    area.      From    6-11 
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February    the    300    cab    jet    max      shifted    northward   away   from   the 


aat a    area. 


TABLE    II 
Calculated   Subsidence    Ratss 


RUN 

DATE 


6-7         6- 7         7-3         s-9 


5  6  7  8  9 

1-9       11-12     14-15    21-22    22-23 


METHOD 

H            -.30      -.70      -.65       -.60       -.55      -.55       -.75  -.31       -.50 

Q                  ♦      -1.72            +         -.39      -.06    -1.17       -.55  -.41       -.33 

A            -.54       -.46    -1.10             +               ♦         -.37       -.33  +          -.50 

*   subsidence    rates    in   :m/5. 

**    H   hind  cast,    Q   specific   humidity,  A    adiabatic 


From  9-12  February,  subsidence  rates  remained  at  -.0055 
m/s .  Also,  no  change  in  strength  or  position  of  the  250  mb 
jet    occurred    between    11    a  .id   12   February. 

A  subsidence  rate  of  -.0075  m/s  was  letermined  for  the 
period  between  14  and  15  February  which  was  an  increase 
since  the  12th.  The  300  mb  jet  strengthened  and  shifted  to 
the  south  toward  the  data  area  from  12-15  February,  and  the 
250    mb   jet    increased    in    strength    from    14-15    February. 

During  the  21-23  February  period,  subsidence  rates 
increased    from   -.003  1      to    -.3050    m/s.         During      this    period, 
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both  the  300  mb  and  the  250  mb  jet  saxina  moved  closer  to 
the  data  area  with  the  30  0  mb  jet  st: inqi hening  and  the  250 
mb  jet  weakening. 
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Figure   5 


Hadiey   Cell  and   Subtropical    Jet    Position   Shewing 
Ar^a   of   Maximum   Downward   Motion,    after    Palmer,    and 
Newton     (1  969)  . 


IV.      THE   MODEL   PERFORMANCE 

A.        EUCT    PREDICTION 

Each  sounding  was  run  with  I  REPS  to  establish  a  baseline 
to  assess  the  model's  ability  to  predict  changes  in  ducting 
conditions.  It  was  then  possible  to  compare  the  model  pre- 
dicted duct  heights  with  thoss  observe?..  Seven  different 
pairs  of  soundings  were  examined,  with  =  ach  pair  separated 
by  24  hours.  Two  of  the  24  hoar  periods  had  2  simultaneous, 
spatially  separated  soundings  at  the  end  of  the  periods.  On 
one  occaisicn,  one  of  the  second  soundings  was  taken  at  a 
location  4  25  nautical  miles  to  ths  southeast  of  the  first 
sounding . 

The  observed  winds  at  the  start  and  and  of  the  period 
were  usad  and  interpolated  linearly  within  the  period  to 
specify   the      wind    for   the      prediction    period.  Sea    surface 

temperature  at  the  latitude  and  longitude  of  the  initial 
sounding  and  at  the  latitude  and  longitude  of  the  verifying 
sounding    were   also    in-.erpclated    within    ths    period. 

Subsidence  velocities  weca  determined  by  hindcasting 
such   that      the   predicted    inversion      height    agreed      with   that 
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observed  in  the  verifying  sounding.  Predicted  ducts  using 
this  "hindcast"  subsidence  are  compared  with  IREPS  observed 
ducts     in   Table     III    .  The   RMS      errors    of      the   model      are 

approximately  one  half  the  errors  of  persistence.  Persis- 
tence is  considered  to  be  a  "prediction"  of  no  change.  It 
is  frequently  used  as  a  baseline  against  which  predictions 
are  compared  to  determine  if  the  predictive  method  has 
merit.  These  results  indicate  that  when  subsidence  can  be 
accurately  estimated,  the  model's  performance  is  clearly 
superior    to   persistence    in   oredic-ing    due1:    bases    and   -ops. 


TABLE    III 

Comparison    of    Observed    Duct    Heights   With    Model   Predicnei 
Heights   Using  Hindcast    Subsidence    Values. 


RUN 

INITIAL 

VERIFYING 

MODEL 

NUMBER 

SOUNDING 

SOUNDING 

PREDICTION 

1 

788-983 

4  70-877 

345-740 

2 

788-983 

533-988 

500-960 

3 

470-  87-7 

46- 504 

0-370 

U 

47-504 

0-135 

0-170 

5 

4  7-504 

0-304 

0-156 

6 

no    duct 

5-  584 

11 5-375 

7 

7  23-1031 

and      0-112 

248-569 

0-495 

8 

206-503 

no   d uc~ 

50-420 

9 

no    duct 

0-258 

0-250 

TOP 

173 

3ASE 

1  11 

THICKNESS 

185 

RMS    ERRORS:  PREDICTION         PERSISTENCE 

364 
214 
285 

*  heights    in   meters. 
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Figure   6,       Mcdel    Predicted   Ducts     (dashed)     vs.    Observed    Duc~s 
(solid)    for  the    Period   6-8    February    1930. 

For    ?    true   predictive   application,    it    least    3    sequential 
soundings      are    required.  The      first      and    second      sounding 

enable  the  e-^t  imation  of  subsidence.  Ths  model  is  then  ini- 
tiated from  -he  second  sounding  using  the  subsidence 
obtained  from  the  first  and  second.  rha  -bird  sounding  is 
used  to  verify  the  prediction.  Gisason  (1982)  found  "--hat  a 
m=thcd  which  is  based  on  the  w=ll  mixed  specific  humidity  is 
an  accurate  method  to  calculate  subsidence.  However,  assum- 
ing   persisterce    in    subsidence    rates   appears    to    be   ths    best. 
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A  .steady  decrease  in  subsidence  rates  (Table  IT)  during 
the  firs-  week  of  the  data  period  coincided  with  a  steady 
decrease  in  surface  pressure.  The  upper  level  jet  stream 
shifted  away  from  the  data  area  and  weakened  during  this 
period.  This  suggests  that  judicious  modification  of  per- 
sistence-obtained subsidence  may  be  appropriate  in  some 
cases      especially   when  a    trend   has   been   observed. 

From  three  p=riods  having  3  sequential  soundings,  Tables 
IV  and  V  shew  the  results  of  predictions.  In  the  first 
period  (Run  number  3)  ,  an  elevated  layer  became  nearly  sur- 
face oased.  The  model  prediction  was  for  a  surface  based 
duct.  It.  the  second  period  (Run  number  U)  ,  a  low  elevated 
duct  became  surface  based.  The  model  prediction  matched  the 
observation  nearly  perfectly.  In  the  third  case  (Run  number 
9)  ,  a  surface  based  duct  formed  when  -her?  was  no  duct  ini- 
tially. In  this  case,  the  model  predicted  the  formation  of 
an    elevated      due-.  While   the      predicted    due4:    in      the    third 

case  differed  subst  a  rtially  in  both  height  and  thickness 
from  the  observed  duct,  the  model  prediction  was  significant 
because  the  occurance  of  a  duct  was  predicted  from  an  ini- 
tial   non -ducting   condition. 
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For  this  particular  comparis  or. ,  tha  absolute  errors  are 
more  meaningful  than  rms  errors.  Two  of  three  model  pre- 
dictions are  essentially  on  the  mark  while  all  of  the  error 
is  concert  re  ted  in  the  third.  In  the  persistence  cases,  the 
errors  are    fairly   evenly    distributed. 


TABLE    IV 

Comparison   of    Observed    Ducts   with   Model    Predicted   Ducts 
Assuming  Persistant    Subsidence. 


?  UN 
NUMB  EH 

INITIAL 
SOUNDING 

VERIFYING 

SOUNDING 

MODEL 
PREDICTION 

3 
4 
9 

470-  3*77 
46-504 

r.o    duct 

46-504 
0-  135 

0-258 

3-350 

0-140 

490-820 

*   heights 

in   meters. 
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TABLE    7 
Error    Analysis   of   Model.  Prediction    vs.     Persistence 


RUN 


ABSOLUTE     ERROR 
HODEI         PERSISTENCE 


3    RUN    RMS     ERROR 
MODEL       PERSISTENCE 


TOP 
BASE 

THICKNESS 

TOP 

BASS 
THICKNESS 

TOP 

BASE 

THICKNESS 


154 

46 
108 

5 
0 
5 

562 
490 

72 


373 

4  24 

51 

369 

47 

3  23 

253 

0 
253 


336 

284 
75 


337 
246 
240 


*  heights    in    meters. 

B.        CLOUD    PREDICTION 

The  cloud  format icn  is  predicted  on  the  basis  of  concur- 
rent prediction  of  the  lif-ing  condensat iDn  level  (LSL)  and 
the    inversion    height.  If   the    LCL   is      bslDw   the   inversion, 

clouds  occur;  if  it  is  above,  no  clouds  occur-  When  model 
predictions  showed  do  low  level  clouds  forming,  satellite 
imagery  confirmee  -.here  were  no  low  lsvil  clouds  and  than 
the  sounding  site  either  was  clear  or  shnwed  scaA-e:ed  high 
cirrus  clouds,  The  formation  of  low  level  stratus  clouds 
was  predicted  en  3  occasions  over  the  entire  period.  Stra- 
tus clouds  were  predicted  to  form  at  2100  local  time  11  Feb- 
ruary. The  first  satellite  imagery  available  following  this 
time  was  at  0300  13  February,  and  showed  the  presence  of  low 
level   stratus    clcuds     (Fig.    8)  . 
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Fiaure  7.   Model  Pred ictior.  Showing  LCL  Droppina  Beneath  the 
Inversion  at  2100  11  Feb. 

Clouds  were  predicted  to  form  at  0500  15  February  (25 
hours  into  the  30  hour  prediction).  In  this  case  the  LCL 
dropped  slightly  below  the  inversion,  but  less  sharply  than 
in  the  other  cases.  Satellite  imagery  from  15  February 
shows  no  clouds  (Fig.  10)  . 

Satellite  imagery  shows  tha  presence  of  low  level  stra- 
tus clouds  on  the  2000  pass  1°  February  and  increased 
amounts  on  the  0800  2  0  February  pass,  Although  no  soundings 
wsre  available  to  initialize  the  model  prior  to  this  initial 
formation  of  clouds,  on  the  next  available  model  run  clouds 
ware  pradicted  throughout,  the  period  0730  21  February  to 
0700  22  February.   Satellite  imagery  confirm  the  presence  of 
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Fiaure   3 


NOAA    6    Satellite   IR   Imagery    from    }800    13    February 
Showing   the  Presence    of  Stratus    Clouds. 


both   cirrus      and   stratus      clouds    on    both      the    0800      and    2000 
passes   on    21    February    (Fig.    12). 

Although  the  satellite  imagery  analysis  was  by  neces- 
sity, subjec-ive,  some  objectivity  was  retained  by  confirm- 
ing the  initial  analysis  with  a  secon i  analysis  by  an 
independent   source.  Both  analyses   agreed    on      cloud    types. 

Since   the    model      could  be    shown    to   be    incorrect      only    on   one 
occasion,      the   model's  overall    performance    in    predicting   the 
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Figure  9.      Mcdel   Prediction   0300    14   to    0900    15    February. 


'igure    10.      NOAA    6    Satellite   IR    Imagery    from    15    February. 
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Fiqure   11.       Model    Prediction    073  0    21    to     (3  00    22    February 


Figure    12.      NOAA    6    Satellite    IP.    Imagery    3330    21    February. 
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formation  of  low  level  clouds/fog  was  sight  out  of  nine 
cases  correct   or   89%    correct.. 

C.        EXAMINATION    OF    MOTEL    SENSITIVITY 

Surface   observation    of    sea    surface    temperature    (SST)     and 
wind   can    be   obtained      with   an    accuracy    of    +/-    1 °C      and    +/-    2 

knots  re-spec-  ively .  To  test  the  model  sensitivity  to  these 
possible  measurement  errors,  successive  model  runs  were  made 
with  SST  1 °C  higher  and  lower  than  observed.  The  effect  of 
wind  measurement  error  was  evaluated  using  wind  values  2 
knots  higher  than  observed.  The  resulting  duct  predictions 
were  then  compared  with  ducts  predicted  by  actual  observed 
winds   and    SSTs.  Table    VI   contains    the      tabulated   results. 

Noting  that  F.MS  error  in  duct  specifications  due  to  poten- 
tial SST  measurement  ^rror  a.ce  approximately  double  the  RMS 
error  due  to  potential  vine,  measurement  srrcr,  one  might 
reasonably  conclude  SST  determination  is  the  more  critical 
of  the  two  parameters.  It  is  iiiiperativs  to  recognize  that 
the  win!  has  a  much  wider  range  cf  variability  than  the  SST 
and  accurate  specification  of  both  parameters  are  important 
for    good    model    performance. 


T?y  "5  T  %'       ITT 

Sensitivity   of    Model    Duct    Predictions    t  3    SST    and   wind 

Variations, 


RUN  VERIFYING  SST  SST  WIND  WIND 

NUMBER  DUCTS  +1 °C  -1°C  +2    KTS         -2    KTS 


1 

345-740 

4  0  0-800 

3  15-715 

390-8  10 

330-725 

2 

500-960 

780-1190 

4  20-370 

62  5-1080 

545-975 

3 

0-370 

0-395 

0-360 

0-400 

0-365 

U 

0-170 

0-273 

0-1  10 

0-130 

0-140 

5 

0-153 

0-210 

0-1  15 

0-170 

0-130 

6 

115-375 

17  0-570 

0-8  0 

25-3  10 

0-240 

7 

0-495 

290-430 

0-4  4  0 

0-510 

0-450 

8 

50-420 

0-360 

0-3  10 

7  0-450 

0-370 

9 

0-25  0 

0-240 

0-290 

0-2<l0 

0-190 

Combined    RMS    error:       SST  DUCT  WIND 

114  top  55 

104  base  50 

107  thickness        26 

*  heights    in   meters. 

The  model's  demonstrated  sensitivity  to  sea  surface 
temperature  points  to  the  need  for  accurate  3ST  determina- 
tions. Current  methods  to  measure  sen  ::.  irface  temperature 
include:  monitoring  sea  water  injection  temperature  (20-^0 
ft.  beneath  the  surface)  ;  utilizing  27  J.rops  (get  tempera- 
ture at  abou*  1  ft.)  ;  and  the  bucket,  and  thermometer  method 
(get  temperature  of  top  few  inches),  None  of  these  methods 
yield  the  sea  surface  skir.  temperature  which  is  the  quantity 
necessary  tc  determine  the  surface  fluxes.  These  fluxes,  in 
turn,  drive  the  evolution  of  the  mixed  layer,  the  inversion, 
and    thus    ducting   conditions   of   the   lower    atmosphere. 


39 


The  052340Z  February  1980  sounding  vas  used  to  test 
model  sensitivity  to  varying  subsidence  rates*  The  hindcast 
subsidence  was  -.0030  m/s  (Table  II)  which  resulted  in  a 
predicted  3S5  meter  thick  elevated  duct  with  the  top  at  740 
meters  alevation  and  the  base  at  345  meters.  Subsidence  was 
varied  from  this  baseline  valua  and  deviations  in  the  pre- 
dicted duct  height  and  thickness  were  recorded.  This  par- 
ticular sounding  was  chosen  because  the  associated  initial 
and  verifying  ducts  were  elevated,  which  allowed  the  varied 
subsidence  tc  affect  the  position  of  the  predicted  duct  in 
both   the    upward   and   downward   directions. 

Deviations  in  the  duct  predictions  resulting  from  subsi- 
dence specification  errors  of  «■/-  .0030  m/s  are  245  m  a:  the 
top,  26  9  m  at  the  base,  and  25  m  in  thickness  (Table  VII). 
Errors  resulting  from  assuming  persistence  in  duct  heights 
(Table  III)  are  approximately  aqual  no  the  errors  resulting 
from  a  subsidence  specification  error  of  +/-  .0030  m/s. 
E£om  this  sinale  case  comjaari sp_a  ,  it  is  ci^ar  thv:  the  model 
i§  £J2  improvement  over,  persistence  only_  when  sjibjsidence  can 
be    specified    to    withir  +  /-    -0030    m/s. 

The  range  of  subsidence  calculated  iuring  the  short 
period      from    6      to    22      February     1980      in    the      data   area      was 
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-.0030  to  -.0075  m/s  or  -.00525  +/-  .00225  ra/s.  Assuming 
this  to  be  e  normal  range  of  variability  and  the  median 
value,  -.00525  m/s,  tc  be  reprssentative  of  the  climatologi- 
cal    subsidence,      then   the      "climatological"    value   appears    to 

be  less  than  the  +/-.0030  m/s  marginal  utility  criterion. 
This  conclusion,  based  on  a  limited  data  set,  seems  to  indi- 
cate that  a  climatol ogical  subsidence  value  would  provide  a 
duct  height  forecast  which,  on  the  average,  has  some  skill 
ever   persistence   or    a    forecast    of   no    change. 

In  reality,  errors  will  not  ba  confined  to  subsidence 
alone,  but  will  be  caused  by  errors  in  other  inputs  such  as 
wind  and  SST.  While  SST ,  wind,  and  subsidence  errors  may  be 
offsetting  at  times,  there  will  be  times  when  -hey  will  be 
additivs.  When  all  the  errors  are  addi-ive,  the  expected 
duct  height  error  for  SST  plus  wind  measurement  errors  are 
16  9  m  at  the  top  and  154  m  at  the  bottom  (Table  7).  Given 
the  errors  due  to  assuming  persistence  in  duct  heights  of 
364  m  at  the  top  and  21 4  m  at  the  base  of  the  duct  (Table 
IV) ,  the  difference  of  approximately  200  m  at  the  top  and  60 
m  at  the  base  indicate  the  maximum  i  uct  height  error  permis- 
sible due  to  subsidence  specification  error  is,  in  the  mean, 
130    m.         From    Table   VI,         the    subsidence    specification    error 
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expected  to  produce  this  130  m  duct  height  error  is  approxi- 
mately +  /-  .0015  m/s. 

From  the  cases  (Table  III)  in  which  the  best  possible 
estimate  of  subsidence,  the  hindcast  value,  was  used,  the 
PMS  model  prediction  error  averages  143  m  less  than  persis- 
tence error.  It  will  be  shown  later  that  if  subsidence  is 
allowed  to  vary  from  this  "besc  estimate",  a  deviation  of 
about  +/-  .0017  m/s  will  produce  the  additional  143  m  error 
in  height  of  duct  top  and  base  necessary  tc  render  the  pre- 
diction no  better  than  persistence. 

In  conclusion,  evaluations  of  the  model  with  the  data 
indicate  that  with  the  expected  SST  and  wind  measurement 
errors,  subsidence  must  be  specified  to  within  +/-  .0015  m/s 
for  the  model  to  consistently  outperform  persistence.  The 
observed  ranae  of  subsidence  variablity  from  this  limited 
data  set  of  +  /-  .00225  m/s  effectively  rules  out  the  use  of 
a  climatclogical  subsidence  as  a  viable  input  parameter. 
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TABLE   711 
Sensitivity    of  Model    Duct    Predictions    to    Subsidence. 


SUBSIDENCE 

PREDICTED 

DUCT 

PARAMETERS 

(ffl/s) 

TOP 

CHANG 

E         BASE 

CHANG 

rE 

THICKNESS 

CHANG 

.0000 

1020 

2  80 

645 

300 

375 

-20 

-.0005 

970 

2  30 

585 

240 

385 

-10 

-.0010 

910 

170 

535 

190 

375 

-20 

-.0015 

865 

1  25 

480 

135 

385 

-10 

-.00  20 

830 

90 

430 

85 

400 

5 

-.00  25 

785 

45 

380 

35 

405 

10 

-.0030    * 

740 

0 

3^5 

0 

395 

0 

-.0035 

713 

-27 

312 

-33 

401 

6 

-.0040 

675 

-65 

250 

-85 

415 

20 

-.0045 

630 

-1  10 

225 

-120 

405 

10 

-.0050 

600 

-1  40 

180 

-165 

420 

25 

-.0055 

570 

-170 

145 

-200 

425 

30 

-.0060 

535 

-2  05 

110 

-235 

425 

30 

RMS    ERRORS: 

DEVIATION    FROM    REFERENCE  EXPECTED    DEVIATION 

SUBSIDENCE         TOP         BASE  THICKNESS  (Ws    X    24    hours) 

+/-    .0030            245            269  25  259 

+/-    .0025            202            221  22  216 

♦/-    .0020            156             178  23  173 

+/-    .0015            118            128  10  130 

+/-    .0010               79               85  15  86 

V-    .0005               37               34  8  43 

*      reference    values 
**      heights    in    meters 


43 


v»        IACTICAL    APPLICATIONS 

The  first  consideration  in  tactical  application  of  the 
model  is  the  recognition  of  its  capabilities  and  limita- 
tions. It  can  be  used  to  estimate  when  low  clouds  or  fog 
can  be  expected  to  form.  It  can  be  used  to  estimate  the 
expected  position  of  low  elevated  and  surface  based  ducts 
below  1200  m  over  an  18  -  24  hour  period.  However,  it  can- 
not   predict    information   on   upper    level    ducts. 

The  tactical  significance  of  knowing  up  to  30  hours  in 
advance  that  visibility  may  be  reduced  or  a  low  ceiling 
developed  will  depend  on  current  and  planned  operations. 
The  results  presented  in  Chapter  IV  B  indicate  that  the 
model   does    well    in    predicting    low    level   stratus/fog. 

When  fog  or  low  level  stratus  clouds  are  predicted,  the 
following  possible  effects  on  own  and  enemy  forces  in  each 
warfare   area   should    be  considered: 

(1)  maneuverability    restricted 

(2)  lest    or   reduced    visual    signalling    capability 

(3)  reduced   visual   target   detection/identification 

(4)  degraded   optical    equipment    performance 
including   IB    weapons    and    sensors 

(5)  hampered    flight    operations    for   carrier    based 
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aircraft   and   LAMPS 

Tactical  employment  of  Naval  assets  requires  a  knowledge 
of  the  state  of  the  environment  if  thoss  assets  are  to  be 
employed      effectively.  To   that     end,        IREPS      (Integrated 

Refractive  Effects  Prediction  System)  was  developed.  When 
radiosonde  data  are  inpur  to  I? EPS,  a  rafractivity  profile 
is  generated  which  is  utilized  to  define  the  location  cf 
ducts  and  to  assess  radar  coverage,  ESM  and  communication 
ranges. 

Because  scundings  are  normally  taken  only  once  or  twice 
nail}/  from  aircraft  carriers  and  the  atmospheric  boundary 
layer  undergoes  change,  -hers  exists  a  need  to  be  able  _o 
predict  changes  in  t he  refracti vity  profile  between  sound- 
ings. This  prediction  can  be  oritically  important  to  battle 
croups  in  the  positioning  of  both  surfaoe  and  air  assets, 
end    in  the    intelligent   management    cf    -he    2MC0N   plan. 

The  duct  prediction  capability,  as  detailed  in  Chapter 
IV  A,  provides  a  significant  planning  tool.  Tactical  useage 
cf  the  duct  predictions  stems  from  the  effects  of  ducts  on 
EM  propagation.  Prediction  of  a  surface  based  duct  indi- 
cates probable  extended  ranaes  for  transmitter-receiver 
antenna    pairs      in   the      duct   and      possible    holes      or    gaps      in 
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coverage  just  above  the  duct  for  antennas  in,  near,  or  below 
the  duct. 

The  predicted  M  profile  may  be  used  as  an  environmental 
data  set.  to  run  IREPS  and  coverage  diagrams  generated  for 
specific  emitters.  This  procedure  was  used  with  a  hypothet- 
ical surface   search  radar  under  observed   refractive  condi- 


tions a.'.'.I  model   predicted   refractive  conditions, 
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version  '.7  (unclassified)  was  u  sea  :.o  generate  the  coverage 
diagrams.  Radar  parameters  for  the  hypothetical  radar  were: 
Antenna  heigh".  110  ft.,  antenna  type  sin(x)/x,  vertical  beam 
width  10 °r  elevation  angle  0°,  free  spacs  range  55  nautical 
miles,    ana    frequency    5  Ghz. 

Clearly  shown  in  Figs.  14  and  15  is  the  model's  value  in 
predicting  the  occurance  of  low  iavei  3ucts.  The  radar 
coverage  In  Fig.  14  is  significant  in  that  the  initial  pro- 
pagation conditions  were  for  surface  target  detections  out 
to  about  25  nm  and  the  predicted  conditions  were  for  detec- 
tions well  in  excess  of  100  nm.  The  pradicted  propagation 
conditions    wer a    in   close    agreement   with   those   observed. 

An  example  of  how  the  modal  might  be  used  operationally 
is      as  follows: 

(1)    radiosonde    is    launched   and   surface    observations 
of    SST   and   wind    are  recorded. 
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Figure   13.       Initial    Surface   Search   Radar    Coverage   0420    7    Feb 
80. 


Figure    14.       Surface    Search  Radar    Coverage    0333    8   Feb    80 
Based    on    Actual   Sounding    (1)     and   Model 
Prediction    (r)  . 
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Figure    15.       Surface    Search   Radar    Coverage    3240    9    Feb    80 
Eased   on    Actual   Sounding     (1)     and    Model 
Prediction    (r)  .. 

(2)  radiosonde   data    are   input   to    IREPS    and   coverage 
diagrams   are    generated   to   make   current   assessment. s 

(3)  radiosonde   data   and  surface   observations   are   input 
to    the    NPS    model    and   a    30    hour   prediction    is    made. 

(4)  the   predicted    .1    profile    is  input    to    IREPS    and 
predicted  coverage   diagrams   are   generated. 

(5)  operational    plans,    incorporating   forecast    ducting 
effects,    are    proposed. 

(6)  tactical   decisions,    using   the    environment    to    the 
advantage  of    the    force,    are    made. 
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In  a  scenario  similar  to  that  observed  in  the  Arabian 
Sea  7  to  9  February  1980  when  an  elevated  duct  became  sur- 
face based,  niany  tactical  considerations  ire  involved.  Sup- 
pose the  battle  group  wanted  to  avoid  detection  by  an  enemy 
force  in  the  region.  Forewarned  that  surface  based  ducting 
was  exvzct&c.  :o  occur  in  12  to  18  hours  and  that  JHF  radio 
transmissions,  normally  useful  between  ships  only  up  to  25 
to  30  nautical  miles,  could  be  detected  at  ranges  in  excess 
cf  200  nautical  miles  due  to  ducting,  the  battle  group  com- 
mander would  be  able  to  impose  a  more  restrictive  EMCON  con- 
dition in  ad£gua+e  time  to  prevent  hostile  intercept  of  the 
group's  transmissions.  Alternatively,  if  the  battle  group 
commander  expected  to  be  opposed  by  active  surface  search 
radar,  he  could:  initiate  a  preeraptiva  attack  before  his 
presence  was  exposed;  steam  out  of  range;  or,  relax  EMCON  at 
the    appropriate   time. 

For  a  battle  group  in  a  multithreat  environment,  early 
threat  detection  is  critical  to  defeat  the  threat  at  minimum 
cost  to  the  battle  group's  offensive  capability.  The  12  to 
18  hour  advance  notice  of  surface  based  ducting  provides  the 
tactical  planner  the  -ool  necessary  to  use  his  limited  sur- 
veillance   assets      with  the      efficiency   needed      for   sustained 
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operations  and  the  e ffectivenes s  require!  for  early  threat 
detection.  Surface  chips  could  be  positioned  in  such  a  way 
as  to  take  maximum  advantage  of  the  extended  surface  search 
radar  ranges  afforded  by  the  sarface  based  duct.  Airborne 
surveillance  assets,  freed  from  task  of  surface  surveil- 
lance, could  concentrate  their  efforts  toward  detecting  and 
tracking  air  targets  .  The  net  result  is  increased  force 
effectiveness.  Depending  on  the  battle  group's  speed  and 
the  picket  station  assigned,  the  surface  ship  may  reguire 
several  hours  to  reach  the  assigned  station,  thus  making  the 
lead   time    provided   by   the   prediction    all   the    more   important. 

The  timing  of  war  at  =ea  strikes  to  coincide  with  pre- 
dicted occurrence  of  a  surface  based  duct  could  give  consid- 
erable advan-.aae  to  the  strike  aircraft  by  enabling  them  to 
use  predicted  holes  ;'ust  above  the  duct  to  penetrate  enemy 
surveillance  and  to  use  the  duct  for  targeting  and  jamming. 
In  such  a  case,  a  12  no  18  hour  prediction  provides  time  to 
plan  mission  tactics,  arm  aircraft,  and  conduct  mission 
briefs. 

When  sonobuoy  patterns  are  laid  and  monitored  by  LAMPS, 
one  constraint  on  placement  of  the  bouys  is  LAMPS'  on  sta- 
tion   endurance.      When   a    surface    based   duct    occurs,    sonobouys 
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car.      be      monitored      at  creator      distancss.  Prediction      in 

advance  of  surface  based  ducting  provides  the  lead  time  nec- 
essary to  plan  and  lay  bouy  patterns  at  distances  from  the 
battle  group  which  take  advantage  of  the  extended  monitoring 
range. 
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VI.       CONCLUSIONS    AND    RECOMMENDATIONS 

The  Tactical  Environmental  Support  System  (TESS)  ,  cur- 
rently under  development,  will  incorporate  various  satellite 
derived  products,  and  single  station  assessment  systems  such 
as  ISEFS  and  the  NPS  mixed  layer  model.  This  will  place 
srate  of  the  art  environmental  sensing  and  predictive  capa- 
bility at  the  disposal  of  those  who  need  the  information: 
the    operational    commanders. 

Eased  on  evaluations  with  a  limited  data  set,  the  NPS 
model  has  been  demonstrated  to  perform  w?ll  both  in  pre- 
dictions of  lew  level  cloud  forma-ion  and  In  predictions  of 
low      elevated   and      surface   based      ducts.  Further   work.      is 

reguired  to  determine  the  model's  rangs  Df  applicability 
both    geographically      and    seasonally.  It    Ls      believed    that 

when  subsidence  can  be  specified  with  an  accuracy  of 
+/-.0015  m/s,  the  model  can  ba  used  to  accurately  forecast 
low      level      ducts.  Two        additional      recommendations      are 

detailed    in   the    following    paragraphs. 

In  tactical  situations  where  the  outcome  may  depend  on 
accuracy    of   duct    predictions,      a      method    is    needed   to    verify 
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the  prediction.  A  radiosonde  Launch  may  be  impractical  due 
to  EMCON  conditions.  An  acoustic  souniar  can  be  used  to 
measure  the  height  of  the  inversion  which,  as  shown  in  Fig. 
3,  corresponds  to  the  top  of  the  duct.  Davidson  ,  et  al 
(1982),  while  conducting  marina  atmospheric  boundary  layer 
research,  utilized  the  economical  acoustic  sounder  on  R.  V. 
ACANIA  to  monitor  inversion  height.  Using  a  combination  of 
daily  radiosondes,  NPS  model  duct  predictions,  and  periodic 
acoustic  soundings  to  track  the  inversion  height  wcild  pro- 
vide improved  capability  to  use  the  environment.  It  is 
therefore  recommended  that  similar  acoustic  sounders  be 
installed   on    aircraft   carriers. 

Determining  appropriate  temperature  and  humidity  pro- 
files luring  initialization  of  the  model  is  a  possible 
source  of  error.  There  are  two  schools  of  thought  which 
differ  on  the  appropriate  height  for  the  inversion  initiali- 
zation. Sounding  data  typically  show  the  inversion  as  a 
layer  100  to  200  meters  thick,  while  the  model  reguires  a 
single  height  to  be  specified.  One  school  suggests  that  the 
bottom  of  the  inversion  layer  is  best  while  the  second  sug- 
gests that    the    midpoint    of   the    layer    is   bids':    reasonable. 
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The  work  in  this  thesis  used  the  former  method.  The 
model  predictions  (Table  III)  show  a  bias  in  that  nearly  all 
of  the  predicted  duct  heights  are  lower  than  those  observed. 
This  seems  to  indicate  that  inversion  heights  should  be 
digitized  at  the  midpoint  of  the  layer  instead  of  at  the 
bottom.  It  is  recommended  that  further  work  be  dene  in  this 
area   to   establish  a    mere    objective   method    of   digitization. 

The  NPS  atmospheric  boundary  layer  model  could  provide 
the  fleet  a  predictive  capability  and  tactical  planning 
tool.  Consideration  should  b9  given  10  introducing  it  into 
the  fleet  for  operational  evaluation.  The  need  to  under- 
stand the  model  assumptions  in  order  to  assess  model  appli-  , 
cability  in  a  given  environmental  situation,  the  need  for 
accurate  subsidence  estimates,  and  the  requirement  for  sur- 
face wind  predictions  preclude  model  installation  except 
where  gualified  personnel  are  assigned.  In  this  regard,  the 
model  should  be  installed  at  locations  (aircraft  carriers) 
where  there  are  Geophysics  Officers  who  can  interpret  and 
use  the  model  as  a  tool  along  with  other  meteorological 
information. 
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APPENDIX    A 


RADIOSONDE    TIMES     AND    LOCATIONS 


TABLE   VIII 
Radiosonde   Data    Set. 


SOUNDING 

SHIP 

1 

NIMITZ 

2 

NIMITZ 

3 

CORAL    SEA 

4 

NI  MI  TZ 

5 

CORAL    SEA 

6 

NIMITZ 

7 

NIMITZ 

8 

NIMITZ 

9 

NIMITZ 

10 

NIMITZ 

11 

NIMITZ 

12 

CORAL    SEA 

13 

NIMITZ 

14 

CORAL    SEA 

*    d 

ate- time    qr 

DT3 


052340Z 
070020Z 
070030Z 
072333Z 
082240Z 
082333Z 
1  100  14Z 
112346Z 
1323  07Z 
1500  13Z 
210300Z 
212035Z 
212244Z 
2223 15Z 


LATITUDE   LONGITUDE 


23° 
22° 
170 

23° 
230 
930 

23° 
230 

23° 
23° 
24° 
23° 
23° 
23° 


12 

33 
40 
07 
00 
13 
1  7 
10 
03 
1  1 
01 
06 
1  8 
38 


061 
061 
06  5 
061 
06  1 
061 
06  0 
061 
064 
061 
058 
061 
060 
060 


°09 
°36 
°51 
°00 
oil 
°00 

059 

°01 

°30 
°07 
°32 

°58 


9>  E 


qroups   February    1930 


TABLE    IX 
Radiosondes    Used    in  Model    Performance    Runs, 


RUN       NUMBER 

1 

2 

3 

4 

5 

6 

— T 
f 

3 

Q 

SOUNDINGS  : 

INITIAL 

1 

1 

2 

4 

4 

7 

9 

1  1 

12 

VERIFYING 

2 

3 

4 

5 

6 

8 

10 

12 

14 
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